chiral spin-wave edge modes in dipolar magnetic thin films 
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Based on a linearized Landau-Lifshitz equation, we show that two-dimensional periodic allay of 
ferromagnetic particles coupled with magnetic dipole-dipole interactions supports chiral spin-wave 
edge modes, when subjected under the magnetic field applied perpendicular to the plane. The 
mode propagates along a one-dimensional boundary of the system in a unidirectional way and it 
always has a chiral dispersion within a band gap for spin- wave volume modes. Contrary to the well- 
known Damon-Eshbach surface mode, the sense of the rotation depends not only on the direction 
of the field but also on the strength of the field; its chiral direction is generally determined by the 
sum of the so-called Chern integers defined for spin-wave volume modes below the band gap. Using 
simple tight-binding descriptions, we explain how the magnetic dipolar interaction endows spin- wave 
volume modes with non-zero Chern integers and how their values will be changed by the field. 

PACS numbers: 



I. INTRODUCTION 



Spin waves are collective propagations of precessional 
motions of magnetic moments in magnetic materials. 
Magnonics research investigates how the spin wave prop- 
agates in the sub-micrometer length scale and sub- 
nanosecond time scale Especially, the propagation of 
spin waves in periodically nanostructured magnetic ma- 
terials dubbed as magnonic crystals^ — are of one of its 
central interests. Owing to the periodic structurings, 
the spin wave spectrum in magnonic crystal acquires al- 
lowed frequency bands of spin wave modes and forbidden- 
frequency bands dubbed as magnonic band gap.^ Like in 
other solid-state engineering such as electronics, photon- 
ics and plasmonics, the main application direction is to 
explore ability of spin waves to carry and process infor- 
mation. Compared to others, magnonics has a much bet- 
ter prospect for miniaturization of the device, because the 
velocity of a spin wave is typically several orders slower 
than those of light and electrons in solids. 

Recently, the authors proposed a spin-wave analog of 
integer quantum Hall (IQH) state, which has unidirec- 
tional edge modes for spin-wave propagation. IQH state 
is a two-dimensional electron system with broken time- 
reversal symmetry, which supports unidirectional elec- 
tric conducting channels along the boundaries (edges) of 
the systemJ^: The number of the unidirectional (chiral) 
edge modes is determined by a certain kind of topological 
number defined for bulk electronic states, called as the 
first Chern integerJ^— Based on a linearized Landau- 
Lifshitz equation, we have generalized the Chern integer 
well-established in quantum Hall physics into the context 
of the spin wave physics, to argue that non-zero Chern 
integer for spin-wave volume-mode bands results in an 
emergence of chiral spin-wave edge modeJ^ 



The proposed edge mode has a chiral dispersion with a 
band gap for volume- mode bands, which supports a uni- 
directional propagation of spin degree of freedom for a 
frequency within the gap. The sense of rotation and the 
number of the chiral mode is determined by the topologi- 
cal number for volume- mode bands below the gap, which 
itself can be changed by closing the band gap. These fea- 
tures allow us control the chiral edge modes in terms of 
band-gap manipulation, which could realize novel spin- 
tronic devices such as spin current splitter and spin- wave 
logic gates. To have these devices in real experimental 
systems, however, it is quite important to have a number 
of actual magnonic crystals, in which spin- wave volume 
mode bands take various non-zero Chern integers. 

From its electronic analogue,^ -'^^ it is expected that 
non-zerzo Chern integers for spin-wave volume-mode 
bands result from strong spin-orbit coupled interactions, 
such as magnetic dipole-dipole interaction. Namely, hav- 
ing an inner product between spin operator and coor- 
dinate operator, the magnetic dipolar interaction locks 
the relative rotational angle between the spin space and 
orbital space, just in the same way as the relativistic spin- 
orbit interaction does in electronic systems J^iii As a re- 
sult of the spin-orbit locking, the complex- valued charac- 
ter in the spin space (i.e. one of the three Pauli matrices) 
is transferred into wavefunctions in the orbital space. Es- 
pecially, in the presence of finite out-of-plane ferromag- 
netic moments in the spin space, the symmetry argument 
allows the Chern integer for volume-mode bands to have 
non-zero integer- value. In the recent work, employing a 
standard plane- wave theory, we have showed that a two- 
dimensional (x-y) bi-component magnonic crystal under 
an out-of-plane field (along the z-direction) acquires spin- 
wave volume- mode bands with non-zero Chern integers, 
when magnetic dipolar interaction dominates over short- 
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ranged isotropic exchange interaction. From the state-of- 
art nanotechnology, however, it is not easy to synthesize 
the proposed bi-component magnonic crystal experimen- 
tally. Moreover, the proposed model is not simple enough 
to see how magnetic dipolar interaction leads to non-zero 
Chern integers for spin-wave volume-mode bands. 



In the present paper, we introduce much simpler thin- 
film magnetic models, which also support spin- wave vol- 
ume modes with non-zero Chern integers and chiral spin- 
wave edge modes, under the field normal to the two- 
dimensional plane. Based on the models, we show that 
the chiral edge modes have frequency- wavelength disper- 
sions within a band gap for spin- wave volume modes, 
and their chiral directions are determined by a sign of 
the Chern integer for a spin-wave volume mode below 
the gap. Using a simple tight-binding model composed 
of 'atomic orbit als', we further argue that the level in- 
version between the parity-odd atomic orbital (such as 
p-wave orbital) and parity-even atomic orbital (such as 
s-wave orbital) leads to a band inversion, which endows 
spin-wave volume-mode bands with non-zero Chern in- 
tegers. We expect that these findings would give useful 
prototype models for future designing of more realistic 
magnonic crystals which support topological chiral spin- 
wave edge modes. 



The organization of the paper is as follows. In the next 
two sections, we introduced the models studied in this pa- 
per (sec. II and Fig. [1]) and formulate our problem and 
summarize a calculation procedure of spin- wave band dis- 
persions and the topological Chern integers (sec. III). In 
section IV, we show how chiral spin- wave edge modes ap- 
pear and how they change their directions on increasing 
the field. The results shows that the sense of the rotation 
of chiral edge mode is indeed determined by the sign of 
the Chern integer defined for the spin- wave volume mode. 
In section V, we introduce a tight-binding description of 
linearized Landau-Lifshitz equations in the context of the 
present models. In sec.VA, we first clarify spin- wave ex- 
citations within a unit cell in terms of a total angular mo- 
mentum variable. Based on the 'atomic orbitals' thus ob- 
tained, we construct a tight-binding model for a square- 
lattice model (sec. VB) and for a honeycomb-lattice 
model (sec. VC). Using this tight-binding model, we 
explain how a level inversion between different 'atomic- 
orbital' levels leads to an inverted spin-wave band with 
non-zero Chern integers and how the signs of the Chern 
integers are changed as a function of the field. To see how 
the proposed chiral spin-wave edge modes could be seen 
in experiments, we simulate the Landau-Lifshitz- Gilbert 
equation for the square-lattice model near the saturation 
field (sec. VI). The section VII is devoted for summary 
and future open issues, in which we also discuss the ef- 
fects of disorders associated with lattice periodicity and 
shape of the boundaries. 



II. MODEL 

In this paper, we consider two-dimensional periodic ar- 
rays of ferromagnetic islands. We assume that each fer- 
romagnetic island behaves as a single spin and ferromag- 
netic islands are coupled via magnetic dipolar interac- 
tion. In fact, two-dimensional periodic lattice structures 
composed of submicrometer-scale ferromagnetic islands 
have been fabricated experimentally, in which they con- 
firm that each island behaves as a giant single spin under 
some circumstances 

To have volume- mode bands with finite Chern integers, 
we generally need multiple-band degree of freedom within 
a unit cell of magnonic crsytal. To this end, we consider 
two models; decorated square-lattice model and honey- 
comb lattice model (see Fig. [1]). A basic building block 
of both models is a cluster of ferromagnetic islands. For 
the decorated square lattice model, four ferromagnetic 
islands form a circle-shape cluster which encompasses a 
site of the square lattice. For the decorated honeycomb 
lattice model, 3 neighboring islands form a circle which 
encompasses either an A-sublattice site or a B-sublattice 
site of the honeycomb lattice. 

Experimentally speaking, it is also quite likely that 
a submicrometer-scale ferromagnetic island has a num- 
ber of low-energy excitation modes having different spin 
textures within the island. Such modes can be also re- 
garded as multiple-band degree of freedom, so that a sys- 
tem with only one ferromagnetic island within a unit cell 
of magnonic crystal^^ could also have a chance to provide 
volume mode bands with finite Chern integer and asso- 
ciated chiral edge modes. We expect that the theoretical 
results obtained in the present model study would also 
provide useful starting points for further studies on such 
systems 

III. FORMULATION 

For the models introduced above, we first determine a 
classical spin configuration which minimizes the following 
magnetostatic energy; 

E = -^ {AVf Ma{ri)fab{ri - rj)Mb{rj) 

+ HAVj2M,{rj). (1) 

i 

where specifies a spatial location of a ferromagnetic 
island (classical spin). For simplicity, the norm of each 
spin is fixed; \M{rj)\ = M^. The magnetic dipole-dipole 
interaction is given by a 3 by 3 matrix, 

with a, 6 = x^y^z. The summation over z,j in eq. ([T]) 
are taken over all ferromagnetic islands, while the sum- 
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FIG. 1: (Color online) Periodic array of ferromagnetic is- 
lands (gray circles) (a) decorated square-lattice model, (b) 
decorated honeycomb-lattice model. Each gray point stands 
for a ferromagnetic island (volume element is AV), which we 
assume to behaves as a single big spin whose moment is fully 
saturated (Ms). The spins are coupled via magnetic dipole- 
dipole interaction. We took |ecc| = \ey\ = 2.4, 2r = 1.2, 
Ay = 1.70 and Ms = 1.0 for the square-lattice case, while 
|ei| = lesi = leal = 2.4, 2r = 1.2, AV = 1.0 and Ms = 1.0 
for the honeycomb-lattice case. The primitive translational 
vector ttyu {/J. = x^y) are defined as (ix — 

and ay = Cy for 

the square-lattice model and as ax = 63 — 62 and ay — 61—62 
for the honeycomb-lattice case (see text). 



FIG. 2: (Color online) Top- view of classical spin configu- 
rations for the decorated square-lattice model. The field is 
lower than the saturation field, so that spins have finite in- 
plane components, forming a vortex structure. 



A. classical spin configuration 

1. square-lattice model 

For the decorated square lattice case, we found that 
every four spins within a circle-shape cluster form a same 
vortex, 

Mo(r = r(c0.,50j) = Ms{-s^se^, s^ce^.c^), (4) 



mation over a^b = x^y, z were omitted. A corresponding 
Landau-Lifshitz equation reads 

dtMaivi) = eabc{H6b^z 

- AV ^/6d(r, - rj)Md{rj))M,{r,). (3) 



with Oj = (j = I,-- - ,4) and (se^ce) = (sin^, cos6>), 
such that the classical spin configuration Mo{r) re- 
spects the translational symmetries of the square lattice, 
Mo(r + a^) = Mo(r) (see Fig. [2j). A finite out-of-plane 
component {(f ^) is induced by the field. Above the 
saturation field {H > He = 1.71), all the spins become 
fully polarized along the field {(p = tt). 



AV is a volume element for each ferromagnetic island. 
From dimensional analysis, one can see that a satura- 
tion field and resonance frequency of spin-wave excita- 
tions are scaled by MgAV/l^^ where / is a characteristic 
length scale for the periodic structuring within the two- 
dimensional plane, e.g. radius (r) of the circle-shape M- 
spin cluster. In the following, we take this value to be 
around 1; Mg = 1, AV = 1.0, 2r = 1.2 for the square- 
lattice case and Ms = 1, AV = 1.7, 2r = 1.2 for the 
honeycomb-lattice case. 



2. honeycomb-lattice model 

For the decorated honeycomb lattice case, the classical 
spin configuration below a saturation field {H < He = 
0.57) breaks the translational symmetries of the lattice, 
while that above the field is a fully polarized state re- 
specting the translational symmetries of the honeycomb 
lattice. For simplicity, we only consider spin- wave exci- 
tations above the saturation field for the decorated hon- 
eycomb lattice case. 
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B. linearized Landau-Lifshitz equation 



with an Hermite matrix H, 



Starting from the classical spin configurations thus 
obtained, M'o(r), the Landau-Lifshitz equation is lin- 
earized with respect to a small transverse fluctuation fleld 
m^(r), with M{r) Mo{r) + m±{r) and m±_ ± Mq. 
In terms of a rotated frame with a 3 by 3 rotational ma- 
trix R{r), with which M'o(r) is always pointing along 
the 2;-direction, R{r)Mo{r) = MgCz and R{r)m±{r) = 
m{r)^ the linearized equation of motion for the trans- 
verse moments takes the form; 

+ MsAV e^^ Ux{ri, rj)mx{rj) (5) 

where m = {mx^rriy^O) and the summation over the 
repeated indices //, z/, A are taken only over y with 
^xy = —^yx = 1- The flrst term in the right hand side 
includes a demagnetization fleld and the external fleld; 



E. 



a{r,)Mo{n) = -AV > J f{r, - r,)Mo(r,) + He, 



where, provided that M'o(rj) gives a local minimum for 
the magnetostatic energy Eq. ([1]), the equality always 
holds true for a certain scalar function a{ri). The dipole- 
dipole interaction in the second term of eq. (|5j) is given 
in the rotated frame; 

f{n,rj)=Riri)firi-rj)R\rj). 

In terms of m± = rrix ± irriy, which are magnon cre- 
ation/annihilation flelds respectively, the equation of mo- 
tion reduces to a following form; 



-idtCFs 



m-{ri) 



a{ri)Ms 



m^{ri) 
m-{ri) 



+ M,AVy( {++;^^'^^-J Wm4r,) 

(6) 

where a 2 by 2 diagonal Pauli matrix as takes +1 for the 
creation fleld (particle space), while take —1 for the an- 
nihilation fleld (hole space). A Green function /a/3(r,r') 
(a, /3 = ±) in the second term takes a form of a certain 
Hermite matrix in the particle-hole space; 

f++{r,r') f+-{r,r') 
/_+(r,r') /__(r,r') 

1 / I « \ ( fccAr,r') hy{r,r') W 1 1 

2 V 1 "V V fy^^^^') fvyir,r') ) \ -i i 

with flpir.r') = fpa{r',r). Accordingly, the prob- 
lem reduces to solving a following generalized eigenvalue 
problem; 



^\ y.,r, ^ m-{rj) 



(8) 

The sum of j is taken over all spins in the systems. Using 
the Cholesky decomposition^ the Hermite matrix can 
be diagonalized by a paraunitary transformation matrix 



HT = asTE 



(9) 



with a proper normalization condition T^a^T = and 
a diagonal matrix E. 

Now that the saddle point solution respects the trans- 
lational symmetries, ^0(7" + a^) = M'o(r), so does 
the Green function and the demagnetization coefficient, 
f{r + a^,r') = /(r,r' - a^) and a{r + a^) = a{r) 
with the primitive translational vectors {ji = x^y). 
Moreover, the classical spin conflguration eq. (j4]) is in- 
variant under the simultaneous C4 rotations in the spin 
space and the lattice space (around z-axis), so that the 
demagnetization coefficient within a unit cell has no spa- 
tial dependence, a{rj) = a. This also holds true for the 
honeycomb lattice case considered. 

With the Born-von Karman boundary condition, the 
eigenvalue problem reduces to a diagonalization of fol- 
lowing Bogoliubov-de Gennes type Hamiltonian for every 
crystal momentum k = {kx^ky); 



idtCTs 



with 



u--k{ri) 



Mu 



and 



and 



m±(r + a^) = ^e±^'^"-ii±,;,(r). 



ihb 



The summation with respect to j (or Vj) in the right 
hand side is taken over a unit cell. For decorated square 
and honeycomb lattice, Mu = 4 and 6 respectively. The 
summation over the translation vectors b are taken over 
sufficiently many unit cells in actual numerical calcula- 
tions, b = nax + may with —50 < n,m < 50. In terms 
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of the Cholesky decomposition, the 2Mu x 2Mu BdG 
Hamiltonian is diagonahzed 



(11) 



with the normahzation condition, ((il)j\(T^\il)j) = (—1)^^ 
where cFj = for particle bands, j = 1, • • • and 
(jj = 1 for hole bands j = Mu + , 2M[/. Provided 
that the spin- wave Hamiltonian is derived from an energy 
minimum of the magnetostatic energy Eq. ([1]), it is guar- 
anteed that eigenvalues for particle bands (j = 1, • • • , M) 
are positive definite Ej^^ > for any while those for 
the hole bands ( j = Mu + , 2Mu) are all negative, 
Ej^k < for any k. In fact, this is true for all the cases 
studied in this paper. 

The eigenvalues in the particle bands, Ej^k {j = 
1, • • • , Mf/), determine wavelength- frequency dispersion 
relations for all the spin- wave volume- mode bands. An 
eigenvector, is a 'Bloch wavefunction' for the corre- 
sponding spin-wave volume-mode band. In terms of the 
Bloch wavefunction, we have calculated the first Chern 
integer defined for each spin- wave band as^ 

Ch, = i{-ir^e^, I d?k{dkJj\cT:,\dk^^j) 
Jbz 





BZ {^jM - Em.kY 



0.2 



j|ili!Hill[fl!i"'"'''"''^^'jj|J'J™^''TOJIUJil{ 

i!l!!I!l!JR^"-^ili 





dkoj 



(12) m(x) x('r) M(X) 




Contrary to the Chern integer defined for electron's wave- 
function^ eq. ([12]) contains the diagonal Pauli matrix 
0-3 between bra-state and ket-state, which takes +1 in 
the particle space while —1 in the hole space. This 
additional structure comes from the fact that magnon 
obeys the boson statistics, which enforces the respec- 
tive BdG Hamiltonian such as eqs. (|8|10p to be diag- 
onahzed in terms of a paraunitary matrix instead of a 
unitary matrix. Due to this paraunitary character in the 
particle-hole space, we can also argue that the sum of the 
magnonic Chern integer over all particle bands always re- 
duce to zero, Yl^^i = which leads to the absence 
of gapless topological chiral spin- wave edge modeJ^ 

In the next section, we have calculated spin- wave ex- 
citations with the open boundary condition along one 
direction (^/-direction) while the periodic boundary con- 
dition along the other (x-dir.); the frequency- wavelength 
dispersions for the spin- wave edge modes are obtained as 
a function of (surface) momentum along the x-direction, 
kx. The dispersions thus obtained alllow us to see the 
propagation direction of the chiral spin- wave edge mode. 
With changing the strength of the field , we have calcu- 
lated spin-wave band dispersions for both volume modes 
and edge modes and the Chern integer for all the volume 
modes. 



FIG. 3: (Color online) (a-d) Side-view of wavelength- 
frequency dispersions of four spin-wave volume-mode bands 
in the square lattice model under fields (H) normal to the 2-d 
plane; (a) H = 0.0, (b) H = 0A7Hc, (c) H = 0.76Hc (d) 
H = O.S2Hc, (e) H = IMHc, (f) H = l.lHc, (g) H = lAHc 
(h) H = 2.35Hc where He = 1.71. The Chern integer for 
red/blue-colored spin- wave bands is —1/ + 1, while other- 
wise. 



IV. RESULTS 
A. square-lattice model 

Results for the square-lattice model are summarized 
in Figs. I3|4I and [H Without the field, the system has 
no magnetization perpendicular to the plane, so that 
the spin-wave Hamiltonian respects both time-reversal 
symmetry, H-k = H^^ and mirror symmetries, e.g. 
H{kx,ky) — '^{k^-ky)' The Chern integer for all the 
four bands are required to be zero by these symmetries 
(Fig. [3](a)), and no chiral spin- wave edge modes are ob- 
served (Fig. [3KA)). With the field along the z-direction, 
these symmetries are gone. 

On increasing the field, there appear a sequence of 
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FIG. 4: (Color online) (A-D) Wavelength- frequency disper- 
sions calculated with open boundary condition along one di- 
rection (^/-direction) and periodic boundary condition along 
the other (decorated square- lattice model); (A) H = 0.0, 
(B,B-1) H = 0A7Hc, (C,C-l,C-2) H = 0.76Hc (D,D-1) 
H = 0.82Hc with He = 1.71. The system along y-direction 
includes 40 unit cell (L — 40). More than 80% of eigenwave- 
functions for red-colored points are localized from y = L — 3 to 
y — while those for blue-colored points are localized from 
y — 1 to y — A. Compared with Fig. [Slja-d), spectra are com- 
prised also edge-mode bands, whose chiral dispersions runs 
across band gaps for spin wave volume modes. 



band touchings between the lowest spin-wave band and 
the second lowest one at the F-point {H = 0.24i^c) and 
two inequivalent X-points {H = 0.67i^c)- As a result of 
these band touchings, the Chern integers for the low- 
est band and the second lowest one become +1 and 
-1 respectively for 0.24 < H/Hc < 0.67 (Fig. [S^b)), 
-1 and +1 respectively for 0.67 < H/Hc (Fig. [3Kc,d)). 
Correspondingly, there appears a chiral spin-wave edge 
mode propagating in the clockwise direction for 0.24 < 
H/Hc < 0.67, whose dispersion runs across a band 
gap between these two spin-wave volume-mode bands 
(Fig. [31(B,B-1)). When the band gap closes and reopens 



FIG. 5: (Color online) (A-D) Wavelength- frequency disper- 
sions calculated with open boundary condition along one di- 
rection (^/-direction) and periodic boundary condition along 
the other (decorated square-lattice model). (A,A-1) H — 
l.Olifc, (B,B-1) H = l.lHc, (C,C-1) H = lAHc (D,D-1) 
H = 2.35Hc with He = 1.71. 



at H/Hc = 0.67, the chiral spin- wave edge mode changes 
its direction from clockwise to anticlockwise (Fig. |4l^C,C- 
2)). The anticlockwise edge mode remains for a relatively 
larger range of the field strength, 0.67 < H/Hc < 1.4 
(Fig. lli:D),Fig. [5]:A),(B),(A-1),(B-1)). 

There is also another sequence of band touchings be- 
tween the third lowest spin-wave band and the highest 
one. They appear at = 0.71i^c (M-point), H = 
0.79 He (two inequivalent X-points) and H = 0.85 He (F- 
point). Correspondingly, the first Chern integers for the 
third lowest band and the highest band become —1 and 
+1 for 0.71 < H/Hc < 0.79 (Fig. [S^c)), +1 and -1 for 
0.79 < H/Hc < 0.85 (Fig.lSjd)), while otherwise. They 
lead to a chiral spin-wave edge mode with anticlockwise 
propagation (0.71 < H/Hc < 0.79; Fig. lli;C,C-l)) and 
that with clockwise propagation (0.79 < H/Hc < 0.85; 
Fig. |31^D,D-1)) between these two volume-mode bands. 

In the limit of strong field, the system becomes effec- 
tively time-reversal symmetric, = H-k (consult also 
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FIG. 6: (Color online) (A-B) Wavelength- frequency disper- 
sions calculated with open boundary condition along one di- 
rection (y-diT.) with the zigzag boundary (decorated honey- 
comb lattice model). The system along the ^-direction in- 
cludes 30 unit cell (L = 30). Eigenwavefunctions for red- 
colored points are localized from y — L — 2ioy — L (> 80%), 
while those for blue-colored points are localized from y — 1 
to y = 3. (A,A-1) H = imHc, (B,B-1) H = 3.5ifc, where 
He = 0.57. 



a perturbation analysis presented in sec. IVC), where the 
Chern integers for all the four spin-wave volume-mode 
bands reduce to zero and the system does not support 
any chiral edge spin-wave which crosses band gaps for 
spin- wave volume- mode bands. Yet there still exist spin 
wave edge modes, which have parabolic dispersions at 
their lowest (or highest) frequency levels and thus con- 
sist of both right-moving mode and left-moving modes 
on a same side of the boundary. (Fig. [5jD),(D-l)). 



B. honeycomb-lattice model 

Results for the decorated honeycomb lattice are shown 
in Fig. [51 Above the saturation field H > He = 0.57, 
the lowest spin-wave volume-mode band and the second 
lowest one are always separated by a finite band gap. 
The Chern integers for these two bands are quantized to 
— 1 and +1 respectively for H > He = 0.57, and a chiral 
spin- wave edge mode with the anticlockwise propagation 
cross the band gap between these two (Fig. [6^A,A-1)). 
On increasing the field, the band gap becomes smaller 
but always remains finite (Fig. [6fB,B-l)). Only in the 
strong field limit, the gap closes and the lowest two bands 
form two massless Dirac-cone spectra at two inequivalent 
iC-points, k = K and with K - Ci = —K • 62 = 
-K' . 63 = • ei = ^ and K • 63 = • 62 = 0, where 
the Chern integers for these two bands reduce to zero 
(see also sec. IVC). In other words, the band gap and the 
chiral spin-wave edge mode which crosses over the gap 



persist even in the presence of large (but finite) field for 
the decorated honeycomb lattice model. 



V. TIGHT-BINDING DESCRIPTIONS 

Although we are dealing with a classical spin prob- 
lem, the calculation procedure so far is totally in par- 
allel with what standard Holstein-Primakoff approxima- 
tioni^Ti^l does in localized quantum spin models based on 
large S expansion (where S denotes the size of a localized 
quantum spin). Finding a classical spin configuration of 
a given localized spin model (on the order of S'^; treating 
spin as a classical spin) corresponds to the minimization 
of the magnetostatic energy, eq. ([T]) (sec. IIIA). Reduc- 
ing a spin model Hamiltonian into a quadratic form in 
terms of Holstein-Primakoff boson field (on the order of 
S) corresponds to linearizing the Landau-Lifshitz equa- 
tion, eq. ([3]), around the classical spin configuration (see 
sec. IIIB). In fact, we diagonalized a quadratic form of 
the spin- wave Hamiltonian, eq. (j8j), to obtain frequency 
levels of spin- wave modes (sec. IV). A tight-binding (TB) 
description introduced in this section is one approximate 
way of doing this diagonalization, which in fact gives use- 
ful physical pictures for results obtained in the previous 
section. 

To construct a TB description for Eqs. (|6|7p . let us 
first decompose the Hamiltonian defined by eq. (|8|) into 
a diagonal part and off-diagonal part with respect to the 
M-spin cluster index; 

with 

(ii'o)ri,r, = -Msadn^vj 

(13) 
(14) 

where r]r,^r, = 1 - ^r,,r,-, V[ri],[rj] = 1 " ^[r,],[r,-] and [Vi] 
specifies a M-spin cluster in which the spin site Vi are in- 
cluded; the 4-spin cluster for the decorated square lattice 
case (M = 4) and the 3-spin clusters for the honeycomb 
lattice case (M = 3); see Fig. H 

Suppose that Hq is diagonalized in terms of appropri- 
ate 'atomic orbitals' localized at each M-spins cluster; 

HqTq = asToasHo. 

Hq is a diagonal matrix, whose elements correspond to 
respective 'atomic-orbital' levels; The column of Tq are 
specified by atomic-orbital index (j, m), cluster index 
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{n^n') and particle-hole index {(Jl^v). The orthonormal 
condition for the new basis is given as TqCTsTq = a 2,. 
In terms of the basis, the original eigenvalue problem, 
Eq. (pT]) . reduces to 

(^0 + ^1)^ = 0-3^:^, (15) 

with T = TqS. The row of S and the row and the column 
of Hi are specified by atomic orbital (j^m), cluster (n,n') 
and particle-hole (/i, u) indices; 

(^1).. , , . = (T^HiTo),. , , (16) 

(^o)(i,n,/i|m,n',^) = ^j,mSn,n' S^^^^Eoj , (17) 



where 'atomic-orbital' levels Eqj being positive definite, 
Eqj > 0. By its construction. Hi has a finite matrix 
element only between atomic orbitals localized at differ- 
ent clusters, which thus stands for inter- or intra-orbital 
hopping terms between clusters. In terms of the cre- 
ation/annihilation fields for the j-th atomic orbital lo- 
calized at the n-th cluster, 7j^/7j,n7 Eq. (p!5|) takes the 
form. 



m n' 
m n' 

\ 



{j,n\m,n') im,n'J ij,n 



(18) 
(19) 



\ 



M=4 



\ 



M=3 



M=6 



FIG. 7: (Color online) A single M-spin cluster for M = 
4,3,6. M spins (gray circle) align along a circle (radius r), 
taking an equal distance between their nearest neighbor spins. 
Each spin has a saturation magnetization (Ms) and a finite 
volume element (AV^). 



with 



{j,n\'m,n') 



(20) 



A. atomic orbitals 

To gain a useful insight on how 'atomic-orbital' levels 
for Hq behave as a function of the out-of-plane field, let 
us first calculate eigenmodes for a single M-spins cluster 
formed by M spins; M spins align along a circle such 
that any neighboring two spins are separated by a same 
distance (Fig. [7j). As an energy minimum of the magne- 
tostatic energy, the M spins form a vortex structure with 
a finite out-of-plane component, 

M{rj =r(c^.,50.)) (-5^s^,,s^%.,c^) (21) 



with Oi 



M y-J 



1, • • • , M) and (so^ce) = (sin 6>, cos 0). 



The saturation field is given as He /Ms = 6Ao(0) -2Ai (0) 
with 



Ao{0) 



1 , _ 1 



647rr" 



647rr^ ^ — ^ se 



AV denotes a volume element of each ferromagnetic is- 
land (spin) and r is a radius of the circle. For H < 
if = Cos-^[-H/Hc] and a = -4Ao(0) + 2Ai(0), while 
(f = TT Siud a = -H/Ms + 2Ao(0) for H > He. Armed 
with these values, excitation modes are obtained by di- 
agonalizing eq. (p!3|) with rj = r{co-^sOj) and Oj = ^ 
(j = 1, • • • , M). With a proper choice of the U{1) gauge 
for m±; 



1 



eq. (p!3|) can be readily diagonalized in terms of the total 
angular momentum nj; 



Hn 



m^{nj) 
m-{-nj) 



0-3 



with 



m±{nj) 



1 ^ 



m_{-nj) J ' 



where nj is defined modulo M {nj = 0, 1, • • • , M 
Hnj takes the form of 

^++(^j) 9+-{nj) \ 
^-+(^j) 9—{nj) J 



(22) 
-1). 



Hnj = -Msaao - M, 
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H/H, 



0.8 1.0 




FIG. 8: (Color online) (a) 'Atomic-orbital' levels as a func- 
tion of the field in a single 4-spin cluster case. Eq. (|13p is 
diagonalized, where the demagnetization field for each site 
(a) comes from the spins in the same cluster. The saturation 
field at which nj = becomes zero is around He — 1.26. 
The level inversion between nj = 1 and nj = is around 
H/Hc = 0.80, while that between nj = 3 and nj = 2 is 
around H/Hc — 0.83. (a-1) 'Atomic-orbital' levels as a func- 
tion of the field calculated from eq. ([13]) in the decorated 
square lattice model. Eq. ()13p is diagonalized, where the de- 
magnetization field for each site includes not only those from 
the spins in the same cluster but also those from spins in the 
other cluster. The saturation field is around He — 1.71, where 
the level of nj = goes below the others, (b) 'Atomic-orbital' 
levels as a function of the field (single 3-spin cluster case) with 
He 0.32. 



with 

( ] = -2iBo{nj)c,as 

\g-+{nj) 9—{nj) J \ > v 

+ {Ao(nj)(-2 + 34) - Ai(nj)(l + 4)}(7o 

- {3^o(nj)4 + Ai(nj)(l - 



and 

AV 1 

i=i ^ 



with qj = . Eigen- frequency with angular momen- 
tum nj takes the following form in the particle space; 

Snj = MsXnj - 2c^MsBo{nj) (23) 

with 

Xnj = ^[-a + 2Ao(nj) + 2Ai(nj)] 

X ^J[-a^ 2Ao{nj) - 6Ao{nj)cl + 2Ai (nj)c2] . 

Figs. [Ha,b) show how the spin-wave excitations for a 
single cluster with M = 3 and M = 4, eq. ([23]) . behave as 
a function of the field respectively. In either cases, doubly 
degenerate modes at the zero field, nj = 1 and nj = 
M — 1 (p^-wave orbitals respectively in the square-lattice 
case; see Fig. [11]), are split on increasing the field, while 
that with nj = (s-wave orbital) decreases its resonance 
frequency, to reach zero at the saturation field H = He. 
Above the field, the 5- wave atomic-orbital level increases 
again, to form a quasi-degeneracy with the atomic-orbital 
level of nj = M — 2 in the large field limit; 

enj^0=enj=M-2^O{l/H). (24) 

The atomic-orbital levels shown in Figs. [8]^a,b) are 
those for a single 4 (3)-spin cluster, where the demagneti- 
zation field stems only from those spins in the same clus- 
ter. Even when demagnetization fields from surrounding 
clusters are included, which is the case with eq. ([T3]) . 
the field-dependence of the atomic-orbital levels behaves 
qualitatively in the same way as in Figs. [8]^a,b). Namely, 
the decorated square lattice model and honeycomb lattice 
model respects the same 4 (3)-fold rotational symmetry 
as in the single 4 (3)-spins cluster, so that only a value 
of the saturation field and specific expression for demag- 
netization field will be modified, e.g. compare Fig. [Sj^a) 
with Fig. [H; a-1). 

From Figs. [8]^a,a-l), notice that there appear a cou- 
ple of level inversions between different atomic orbitals, 
such as the one between nj = and nj = 1 around 
H = O.SOi^c, and the one between nj = 2 and nj = 3 
around H = 0.83i^c- Now that one of these two atomic 
orbitals is always either + Wy or px — ipy-like orbital 
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M(X) X(r) M(X) ' 

FIG. 9: (Color online) Side-view of spin-wave band disper- 
sions for decorated square-lattice model calculated from the 
tight-binding Hamiltonian, eqs. (|15ll6ll7p where only near- 
est neighbor hopping integrals are included. As for atomic- 
orbital levels, eq. (pT|) . and the respective wavefunction To 
used in eq. {[T6|), we used those for the single 4-cluster. 
(a) H = 0.23ifc, (b) H = 0.66ifc, (c) H = O.Tlifc, (d) 
H = 0.79i7c, (e) H = 0.85ifc, (f) H = lAHc, where He de- 
notes the saturation field for single 4-spin cluster, (He = 1.26; 
see the caption of Fig. [8{a)). At (a),(b) and (f), the lowest 
spin- wave band and the 2nd lowest one form band touchings 
at F-point, X-points, and M-point respectively. At (c),(d) 
and (e), the highest spin- wave band and the 3rd lowest one 
form band touchings at M-point, X-points, and F-point re- 
spectively. 



while the others are parity even, it is expected from 
its electronic analogue^^ that these level inversions en- 
dow the spin- wave volume-mode bands constructed from 
these atomic orbitals with non-zero Chern integers. In 
fact, a similar type of the band inversion plays a vital 
role in the emergence of nontrivial topological phases in 
quantum spin Hall insulators ^^i^^ 





1 
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Fig.9(a) (b) (c) (d) (e) (f) H/H^ 



FIG. 10: Chern integers for four spin- wave bands as a function 
of the field, (a)-(f) depicted in the figure correspond to the 
fields at which band touchings occur as shown in Fig. [9lja-f) 
respectively. Note also that the third lowest spin wave band at 
Fig. [9l^e) is mainly composed of nj = 0, while the 2nd lowest 
and the lowest bands are mainly composed of nj = 2,3. From 
(e) to (f), the former band goes below the latter two until 
H = He, while, for H > He. it increases its frequency again, 
in the same way as the s-wave atomic orbital does in Fig.fSlJa). 
For clarity, we call the latter two as '2nd lowest' and 'lowest', 
even though they are not during (e) - (f). 



B. TB model for the square lattice case 

To clarify how the level inversion between atomic or- 
bitals leads to spin-wave bands with non-zero Chern in- 
tegers, we next construct from Eqs. (|15ll6ll7p a tight- 
binding (TB) model for the decorated square lattice. Hq 
and corresponding Tq in Eqs. ()16|17p are replaced by 
those for the single 4-spin cluster. As for only the 
nearest neighbor hopping integrals are included. Such 
approximations may be justified, because the dipolar in- 
teraction decays as with R being a distance be- 
tween two spins; an amplitude of the next nearest neigh- 
bor hopping is roughly 2>/2, 8, and bVb times larger 
than those of the 2nd, 3rd and 4th nearest neighbor hop- 
ping integrals respectively. In fact, band dispersions for 
spin-wave volume-mode bands obtained from this short- 
ranged TB model show qualitatively good agreements 
with those shown in the previous section (compare Figs. [9] 
with Figs. [3]). For example, all the sequences of the band 
touchings described in sec. Ill are identified near similar 
field strengths, when scaled by the respective saturation 
fields; Fig. [9] (a-f). Correspondingly, the Chern integers 
for all the four spin-wave bands take the same sequence 
of the integer values as found in the previous section; 
Fig. [TOl The comparison also suggests that the non- 
parabolic characters of a certain band dispersion around 
the F-point in Fig. [3] stems from long-range hopping in- 
tegrals in Hi, i.e. long-range character of the magnetic 
dipolar interaction, which is consistent with the similar 
feature of the forward volume modes 

The sequence of band touchings between the highest 
and the third lowest spin wave band results from the 
level inversion between the atomic orbital with nj = 1 
and that with nj = 0; Fig.[8fa), while the other sequence 
between the lowest and 2nd lowest spin- wave band is from 
those with nj = 2 and nj = 3; Fig. [8fa). 

To see this, notice first that the atomic orbitals with 
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nj = 0, 1, 2, 3 takes 5- wave, p_ = — ipy, dx2-y2^ and 
p+ = Pa^+zp^y'Wa^e orbitals respectively; Fig.[TTJ Namely, 
Eq. ([22]) suggests that 'atomic-orbital' wavefunctions for 
nj = 1 and nj = 3 take imaginary values ('i') along the 
^/-direction, while take real values along the x-direction. 
Meanwhile, those for nj = and nj = 2 always take real 
values; nj = takes +1 for x-link and ^-link while nj = 1 
takes +1 and —1 for x-link and y-\mk respectively. As a 
result, the nearest neighbor inter-orbital hopping integral 
between n j = 1 and n j = and that between n j = 2 and 
nj = 3 always become pure imaginary along the y-\mk. 
In fact, using symmetry arguments, one can generally de- 
rive from eqs. (jl8|19p a nearest-neighbor hopping model 
composed of nj = and nj = 1 as; 

^01 = (^07d,n70,n + ei7l,n7l,n) 
n 

+ («007o,n70,n+ae^ + ^ill7l,n7l,n+ae J 

+ (C^^0l70,n7l,n+ae. + h.C.) 

+ ( - ^^^Ol7o,n7l,n+ae^ + h-c), (25) 

n cr=± 

with real valued aoo, <^ii and 601 • We have ignored (or 
integrated out within the 2nd order perturbation) those 
hopping terms which involve nj = 2 and nj = 3 and 
those between the particle space and the hole space. Such 
approximations are justified, when the atomic-orbital 
level of nj = and that of nj = 1 are proximate to 
each other and when it comes to those spin-wave bands 
near these levels. 

The highest and the third lowest spin-wave bands 
around H/Hc = 0.78 ~ 0.82 are mainly composed of 
the atomic orbitals with nj = and nj = 1 (compare 
Fig. IHl^a) with Figs. [9l^c,d,e)) and thus can be approxi- 
mately obtained from Hqi. The Hamiltonian in the mo- 
mentum space takes the for m^^i^^ 

TT _ f £0 ^ 2aoo{ck^ ^ Cky) -2iboi{sk, - isky) 
^oi''^ ^ 2iboi{sk^ + isk^) ei + 2aii{ck^ + c^J )^^' 

with (c/c^ i^k^) = (cos kx, sin kx), which gives momentum- 
frequency dispersions for the highest and the second high- 
est spin- wave bands as 

^o,fc = ^^^^^ + (aoo + o.ii){ck^ + Cky) + (27) 

Si,k = ^^"^ + (^^00 + + ^ky) - (28) 

with 
Afc = 

[eo - ei + 2(aoo - aii)(c/c. + c^J]^ + ^^^oi{4^ + 4J 

(29) 

The atomic-orbital level with nj = and that with 
nj = 1 are inverted around H = 0.80i^c, where £0 — ^1 




nj=2 nj=3 

FIG. 11: (Color online) Shapes of the atomic orbitals in the 
decorated square lattice model, that of nj = is s-wave like, 
while those of nj = 1,2,3 are px — ipy-wave, dx2_y2-wa\e, 
Px + ipy-wave like respectively. 



changes its sign from positive to negative; Fig. [Sfa). 
From their orbital shapes, the nearest neighboring hop- 
ping integral between s-wave orbitals should be always 
positive aoo > 0, while that between p_ -orbitals should 
be negative an < 0, which leads to aoo — ciii > 0. 
These two observations mean that, on increasing the 
field, the two bands given by eqs. (|27l28p first form a 
massless Dirac-cone spectrum at /c = (tt, tt) (M-point) 
when £0—^1 = 4(aoo — <^ii), then two massless Dirac- 
cone spectra at /c = (tt, 0) and k = (0, tt) (X-point) when 
So — ei = 0, and finally one massless Dirac spectrum at 
k = (0,0) when £0 — ^1 = — 4(aoo — <^ii)- The band 
touching at the M-point is nothing but that in Fig.[9fc), 
those at the X-points are those in Fig. [9fd), and that at 
F-point corresponds to that in Fig.[9l^e). In fact, analytic 
evaluations of eq. ([12]) for eq. ([26]) show that the Chern 
integer for the highest (third lowest) spin-wave band be- 
comes +1(— 1) for 4(aoo — <^ii) > Sq — Si > and — 1(+1) 
for > £0 — ^1 > — 4(aoo — <^ii), which is consistent with 
Fig. [101 

Similarly, the other sequence of band touchings formed 
by the lowest and second lowest spin- wave bands is ex- 
plained in terms of the two-band models composed by 
nj = 3 and nj = 2 atomic orbitals; 

TT _ f es ^ ^assick, ^ Cky) -2ibs2{sk, - isk^) V^^x 

^^^'^^ ^ 2ibs2{Sk. + iSky) 62 + 2a22{Ck. + CfeJ /^^^ 

Note that £3—62 changes its sign from positive to negative 
near H 0.83i^c (Fig-[8Ka)), while ass — a22 being always 
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negative. This means that, on increasing i7, the lowest 
and second lowest spin- wave bands around H O.S3Hc 
first form a massless Dirac-cone spectrum at k = (0,0) 
when 63—62 = — 4(a33 — a22), then two massless Dirac 
spectra at A; = (tt, 0) and k = (0, tt) when 63 — 62 = 0, 
and finally one massless Dirac spectrum at /c = (7r,7r) 
when £3 — £2 = 4(a33 — a22). The band touching at the 
F-point is nothing but that in Fig. [9l^a), those at the 
X-points are those in Fig. [Ifb), and that at M-point 
corresponds to that in Fig. [9l^f). Noting that 632 has 
the same sign as 601 (see Fig. [TT]) . one can also see from 
the previous evaluation that the Chern integer for the 
second lowest (lowest) spin- wave band becomes — 1(+1) 
for -4(a33 - a22) > S3 - ^2 > and +1(-1) for > 
63 — 62 > 4(a33 — a22), which is consistent with Fig. [TOl 



C. TB model for the honeycomb lattice case 

In the decorated honeycomb lattice model, we have ob- 
served in sec. Ill a finite band gap between the lowest 
spin- wave band and second lowest one, which are con- 
nected by a dispersion of a chiral edge mode. The gap 
and chiral edge mode persists for a sufficiently large field 
H. Based on a tight-binding model, we will employ a 
perturbation analysis from the large field limit and ar- 
gue that the gap closes at two inequivalent iiT-points only 
in the limiting case {\H\ 00), where both the time- 
reversal symmetry and hexagonal spatial symmetry are 
effectively recovered. More accurately, we will show that 
an effective spin- wave Hamiltonian in the large field limit 
respects these two symmetries within the order of 0(1), 
while it starts to break them from 0{1/H). As a result, 
within the order of 0(1), the lowest and second lowest 
spin-wave band compose massless Dirac spectra at the 
iC-points. Once the 0(l/i^)-corrections are included, 
the time-reversal symmetry is broken and the hexagonal 
symmetry (Cey) reduces to its subgroup symmetry (Ce), 
which leads to a finite band gap at the iiT-points. These 
symmetry breakings also endow the two bands with a 
non-zero Chern integer with opposite signs, which re- 
sults in the emergence of chiral edge mode within the 
band gap. 

The perturbation analysis begins with a tight-binding 
Hamiltonian for the honeycomb lattice model, eq. (p!5]) : 

Ho^Hi= HcTo + XVi + AV2 (31) 

where ctq is a 2 by 2 unit matrix in the particle-hole 
space and both Vi and V2 are on the order of 0{1). For 
a bookkeeping, we put A, which can be set to 1 from the 
outset [those terms with A are 0{1)^ those with A^ are 
0{1/H) and those with A^ are 0{1/H'^); see below]. Vi 
consists of on-cluster 'atomic-orbital' levels and hopping 
terms in the excitonic channel, while V2 consists only of 



those in the Cooper channel; 

^'=[t(-'+) J' 

In the large field limit, all the spin- wave excitations re- 
duce to the ferromagnetic resonance (FMR) with its res- 
onance frequency being H. Once the 0(l)-corrections 
(Vi,V2) are included, the FMR resonance is expected to 
split into a couple of spin-wave bands whose bandwidth 
are at most on the order of unity. To see this situation, 
let us erase those terms in the Cooper channel within 
a given order accuracy in l/H, and derive an effective 
Hamiltonian only for the particle space. 

With a matrix satisfying a 3 pa 3 = p\ the transformed 
Hamiltonian takes the form 

i^eff = e-'^'^^P''^ {H + AVi + AV2) e'^^ 
A^ 

= (1 - i\(T3p(T3 - —Cr3p^(T3 H ) 

X (i7 + AVi + AV2)(1 + i\p - yp' + • • • ) 
= H ^ XVi + AV2 - \i(T3p(T3H + iXHp 

- Y^sP^o-3H - Y^P^ + X^a3pa3Hp 

- i\^cr3pcr3Vi + iX^Vip - i\^cr3pcr3V2 

+ zAV2P + (9(A^) (32) 

We choose p such that all the matrix elements in the 
Cooper channels will cancel each other within the order 
of (9(1); 

V2 icr3pcr3H - iHp. 

Or equivalently. 




where n is for the indices in the hole space and n is 
for those in the particle space; (o-3)n,m = ^n,m and 
(^3)n,m = — ^n,m- With this rotated frame, all the ma- 
trix elements in the Cooper channel are at most on the 
order of 0(l/i^); 

Hef^=H^ XVi - — {CT3PCT3V2 - V2P) 

-iX\cT3pcT3Vi-Vip). (34) 

The last two terms have matrix elements in Cooper chan- 
nels. When we further rotate in the particle- hole space 
such that they will be set off by generated terms, these 
two terms simply result in higher order terms, 0(l/i7^), 
while the remaining terms being kept intact. We thus 
drop them by hand, to keep the first four terms as the 
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effective Hamiltonian. On the substitution of eq. (|33|) 
into eq. we then have the effective Hamiltonian to 

the order of 1/i^ as 



eff /n,m 



^E(^2)„,,(V.)_^. (35) 



The superscript '(2)' in the left hand side denotes that 
the effective Hamiltonian is asymptotically exact within 
the 2nd order in A (or within the first order in 1/H). 
The sum with respect to p is taken only over the hole 
space. From eq. (|35]) . one can ready see that, once the 
0(l)-corrections (Vi) are included, the FMR resonance 
localized at H is split into a couple of spin-wave bands 
whose bandwidth are at most on the order of unity. 

Within the order of 0(1), the effective Hamiltonian de- 
rived so far is invariant under the time-reversal operation 
and hexagonal symmetry operations. To see this, let us 
focus on the first two terms of eq. ([35|) . With the atomic 
orbital index (nj = j, m) and cluster index (n,n') being 
made explicit, they take the following form 



(1) 

ff ^ {j,n\m,n') 



.(+,+) 



In the leading order in the inter/intra-orbital hop- 

ping integral between an orbital with nj = j at the n-th 
cluster and that with n j = m at the n'-th cluster is given 

by 



^(+,+) 

(j,n|m,n') 



1 



(36) 



Oi(') {l^'^ = 1,2,3) in the right hand side specifies a spatial 
location of a ferromagnetic spin within a cluster. Within 
a cluster which encompasses an A-sublattice site or B- 
sublattice site at {xruVn)^ we take Oi = ^ or tt + ^ 
respectively, such that the location of the ferromagnetic 
spin is always given by {xn — rsmOi^i/n + cosOi). R 
denotes a spatial distance between a ferromagnetic spin 
specified by (Oi^n) and that by {Oi'^n')] 



R 



Xn — r sm Ui 
Vn^r cos 6i 



Xn' — r sm f/// 
yn' + r cos Oi^ 



Within the order of 0(1), a complex conjugatation of 
hopping integrals can be readily set off by a sign change 
of the orbital angular momentum ul = nj + 1. Namely, 
the complex conjugate of eq. (|36|) is transformed to itself 
by a proper exchange between = +1 {nj = 0) and 
= 2 = — 1 (mod 3) (nj = 1). This in combination 
with eq. indicates that the effective Hamiltonian up 
to the order of 0(1) is invariant under the following time- 
reversal operation; 



V ett / (j,n|m,n') 



(37) 



with a proper basis change 



Q= 1 



which exchanges ul = +1 ijij = 0) and ul = — 1 
(nj = 1), while keeps ul = {) {nj = 2) intact. It is also 
invariant under three mirror operations in the hexagonal 
symmetry, cr^ i, cr^ 2, cr^; 3, 7r-rotation C2 and ^-rotation 

C3; 



eff / (j,cr^,i(n)|m,cr^,i(n 



\ eff / (j^cr^,2 (n) |m,cr^;,2 (n' 

r(i) 



)(j,C2(n)|m,C2(n' 



/rr(l)\ 

l^eff ;(j,C3(n)|m,C3(n' 



respectively with 



— n (fr'^^h nt 



eff J(j',„|m',„')^»n'm, 



r(i)^ 



,n|m',n') ' 



e 3 



e 3 



O 



e 3 



, p 



^v,v (^ = 1,2,3) denotes the mirror operation with re- 
spect to the plane subtended by e^^ ii' = 1, 2, 3; see Fig. [I]) 
and Bz (a unit vector normal to the plane), C2 is the 
two-fold rotation which exchanges A-sublattice and B- 
sublattice and C3 is the three-fold rotation within the 
plane (see Fig. [T]^b)). 

Owing to the hexagonal symmetry, the lowest two spin 
wave bands obtained from comprise two massless 

Dirac spectra at two inequivalent iiT-points. Once the 
0{1/H) corrections are included, e.g. £0 7^ ^1, the time- 
reversal symmetry is lost and the hexagonal {Cqv) sym- 
metry (C2, C3, ^3"^ (7^,1, (7^,2, Cr^,3, Cr^,! , Cr(i,2 , CFd,^) with 

(jd^v = C2 ■ CFv.v reduces to Cq symmetry (C2, C3, Cg""^). 
As a result, Dirac spectra at the i^-points acquire a fi- 
nite mass, which gives to the lowest two bands non-zero 
Chern integers. 

The Chern integers for the lowest two spin-wave band 
can be evaluated from a nearest neighboring (NN) TB 
model. From the symmetry point of view, the NN TB 
Hamiltonian in the momentum space reads. 



^NNTB,fc 
/ ^0 



e2 



^0 



82 J 
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FIG. 12: (Color online) (a) Side-view of spin-wave band dis- 
persions for the decorated honeycomb lattice model under a 
sufficiently strong field (H = 5.0). Because of a small but 
finite band gap at two K-points, the Chern integers for the 
lowest two bands are —1 and +1 respectively, (b) Side-view 
of spin-wave band dispersions calculated from the effective 
Hamiltonian to the order of 0{1). For the 'atomic-orbital' 
levels, we use those for a single 3-spin cluster. To evaluate the 
nearest neighbor hopping integral within the order of C^(l), we 
use Eq. (|36p . When H — 5.0 is added, the resonance frequen- 
cies of the spin- wave bands in Fig. [I2jb) become comparable 
to those in Fig. [12^ a). 
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The first three columns and rows are for the three atomic 
orbitals encompassing an A-sublattice site, while the lat- 
ter three are for those encompassing a B-sublattice site. 
Sj stands for a level for an atomic orbital with nj = j 
{j = 0,1,2). ao, ai, /3, 7^, 7o and 71 are NN inter or 
intra-orbital (effective) transfer integrals, which can be 
evaluated from eq. ([35]) up to 0{1/H). It is clear from 
Eq. (|36|) that, within the order of 0(1), 70 = 71, Q^o = <^i, 
and So = £1, which makes the lowest two bands form 
massless Dirac spectra at the iC-points; Fig. [I2fb). A 
comparison between Fig. [I2ja) and Fig. [I2fb) suggests 
that the present NN TB Hamiltonian can qualitatively 
well reproduce the band structure of the lowest two bands 
in the large field limit, expect for a non-parabolic band 
structure of the lowest band near the F-point. 

Once finite A7 = 70 — 71, Aa = ao — ai and Ae = 
£0 — Si are included, the exchange between ul = +1 



and riL = —I changes the signs of these terms, so that 
the time-reversal symmetry is broken and the hexagonal 
symmetry reduces to the Cq symmetry. These symmetry 
reductions give a finite mass to the Dirac spectra. The 
mass can be evaluated from 2 by 2 Dirac Hamiltonians for 
the lowest two spin- wave bands, which can be obtained 
via k ■ p perturbation around these iC-points; 

^^i^+P = ^(Ae sin2^-3A7 sin^)^3 

3/ o n • 2 ^\ / \ 

+ 2~^^^^ 2;*^^^^^ 



and 



H 



2x2 

k^K'+p 



1 / \ 

--(^Ae sin^ - - 3A7 sin6>jt73 

^ {jl cos^ ^ - /5 sin^ ^) {-Pxcri - Py(T2), 



with 



tan 6> = 3 



71 +72 

£q — £1 



From these Hamiltonians, the Chern integers for the low- 
est and second lowest spin- wave bands are evaluated to 
be (J and —a respectively with 



/ \ 

sign^Ae sin^ - — 3A7 sin^j. 



A substitution of actual numbers into the parameters in 
the right hand side shows that cr = 1, which is consistent 
with previous numerical evaluations in Sec. HI. The non- 
zero Chern integers for the lowest two spin-wave bands 
results in an edge mode with the anticlockwise propa- 
gation, which has a chiral dispersion between these two 
bands. 



VI. MICROMAGNETIC SIMULATION 

To uphold the existence of proposed chiral spin-wave 
edge mode by a standard method in the field, we perform 
a micromagnetic simulation by solving numerically the 
Landau-Lifshitz-Gilbert equation for the square-lattice 
model. We calculate magnetization dynamics by employ- 
ing the 4th order Runge-Kutta method with a time step 
At = 1 ps. Fig. [13] schematically shows an entire system 
studied in the present micromagnetic simulation. It con- 
sists of 4 ferromagnetic nanograins in the unit cell. Al- 
though we took the size of the ferromagnetic nanograin 
to be 5 X 5 X 5 nm^ as a demonstration, the system is 
scalable; the simulation does not include any short-range 
exchange interactions. The saturation magnetization and 
Gilbert damping coefficient of the ferromagnetic grain are 
1.75 J//im and 1.0 x 10~^ respectively. We regard each 
nanograin as a uniform magnet, to assign single spin de- 
gree of freedom to each nanograin. Different ferromag- 
netic nanograins are coupled with one another through 
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FIG. 13: (Color online) Schematic view of a simulated system 
that comprises ferromagnetic nanograin. Although we took 
the size of the ferromagnetic nanograin to be 5 x 5 x 5 nm^ as 
a demonstration, the present simulation is scale free. In the 
begining of the simulation, we apply a pulse field either at the 
center (marked by a black triangle) or around the boundary 
(marked by a black cross). 



the magnetic dipolar interaction. The simulated system 
(0 < X < L and < r < L; Fig. ^ includes 25 x 25 
unit cells. Without the field, the magnetization of each 
grain lies within the plane due to the dipolar interaction. 
Under a large out-of-plane DC field {H^c > 4700 Oe), 
the magnetization becomes fully polarized along the z- 
direction. We took H = 1.02Hc in the present simula- 
tion. 

In order to excite spin wave modes in a broad fre- 
quency range, we apply a pulsed magnetic field within 
a plane with its pulse time = 1 ps and its amplitude 
Hp = 1.0 X 10~^ Oe. The pulse is applied locally at the 
center and around an edge of the system for the purpose 
of exciting volume modes and edge modes, respectively. 
After calculating a time evolution of the magnetization in 
the system, we take a Fourier transformation of the trans- 
verse moment, m+(X, Y", t) = mx{X^ Y", t) + im^(X, Y", t), 
with respect to time; 

n-l 

s+(X,Y,a;) = ^ m+(X, Y, jAT) exp (27ricjjAT) (38) 

with AT = 50 ps and n = 1024. The frequency power 
spectrum, "^XY \^-\-{^^'^^^)\^ obtained by the pulse at 
the center and that by the pulse at the edge are shown in 
Fig. [TWa) separately. Spatial distributions of spin wave 
excitations, |s+(X, Y, cj) | , for each case with different fre- 
quencies uj are shown in Figs. [lU^b)-(g). From them, one 
can see that the spin- wave volume modes and edge modes 
are selectively excited, depending on whether the initial 
pulse field is applied at the center or at the edge respec- 
tively. In the case of the pulse field at the center, we ob- 
serve two band gaps for volume modes; one from 24GHz 
to 30GHz and the other from 37GHz to 46GHz. In the 
case of the pulse at the edge, we observed spin- wave edge 
modes mainly from 24GHz to 42GHz. 




FIG. 14: (Color online) Fourier power spectra of magneti- 
zation dynamics, (a) Frequency dependence of the intensity 
of spin- wave excitations (see text), (b)-(g) Spatial distribu- 
tion of the intensity, which is obtained by the application of 
the pulse field at center (b,d,f) and the edge (c,e,g); (b,c) 
UJ = lOGHz, (d,e) a; = 29 GHz, and (f,g) a; = 31 GHz. 



A key feature of proposed chiral spin-wave edge mode 
is a unidirectional propagation of spin wave densities, 
which is clarified by its frequency-wavelength dispersion 
relation. To obtain such a dispersion relation, we next 
take a Fourier transformation of the transverse moment 
with respect to both space and time. In order to compare 
the result with Figs. I4|5[ we integrate the amplitude of 
the Fourier component with respect to the Y-component 
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of the momentum; 

ky 

exp (ikxX) exp (ikyY). 

A contour plot of A{kx^oj) as a function of kx and uj gives 
a dispersion relation for spin- wave modes. Figure [TSTa) 
shows the dispersion relation for the case of applying the 
pulse field at the center. It resembles those for spin- wave 
volume modes obtained in the preceding model calcu- 
lation at the same parameter regime (Fig. [Sl^e)). Fig- 
ures [TST b) and (c) show the dispersion relations for the 
case of applying the pulse field at the edge. To clar- 
ify propagation directions of those two spin-wave edge 
modes running along the opposite boundaries of the sys- 
tem, we take the Fourier transformation only over the 
upper (or lower) side of the sample L/2 < Y < L (or 
< Y < L/2); the one for the upper side is shown 
in Figure [TST b). while the one for the lower is in Fig- 
ure [T5fc). Both figures clearly indicate the existence of 
two counter-propagating chiral dispersions, each of which 
runs across any line of u = uq (25GHz < ujq < 35GHz) 
once and only once. The results also suggest an existence 
of another spin- wave edge mode from 35 GHz to 42 GHz, 
which has a quasi-parabolic dispersion. Both of these 
edge modes in combination with volume modes shown in 
Fig. [TST a) are consistent with the dispersion relations for 
spin-wave modes obtained in the preceding model calcu- 
lation at the same parameter regime (Fig. [5fA),(A-l)). 

When the Gilbert damping coefficient becomes larger, 
unidirectional propagations of spin density along the chi- 
ral spin- wave edge mode decay faster. Fig. [16] shows spa- 
tial distributions of Fourier power spectra of magneti- 
zation dynamics, |5+(X, F, cj)|, in the presence of larger 
Gilbert damping term {a = 0.001,0.01), where the ini- 
tial pulse field is applied at the edge (Fig. [T3|) and the 
frequency is chosen within the band gap {uj = 29GHz). 
The results suggest that the coherence length is roughly 
25 unit cell size (500nm) for a = 0.001 and 8 unit cell 
size (160nm) for a = 0.01. 

VII. SUMMARY, DISUCSSIONS AND OPEN 
ISSUES 

In this paper, we introduced two simple magnetic thin- 
film models, in which ferromagnetic nanoislands on pe- 
riodic arrays are coupled with each other via magnetic 
dipolar interaction. Under the field applied perpendicu- 
lar to the two-dimensional plane, spin- wave excitations in 
the systems have a chiral spin-wave edge mode localized 
at the boundaries of the systems, whose dispersion runs 
across a band gap for spin- wave volume modes. The sense 
of the rotation of the chiral edge mode is determined by 
a sign of the Chern integer for a spin- wave volume-mode 
band below the band gap. 




FIG. 15: (Color online) Dispersion relations of spin wave 
modes, (a) Dispersion relation obtained by the application 
of the pulse field at the center. (b),(c) Dispersion relation 
obtained by the pulse at the edge. The Fourier transformation 
is taken only over the upper side {Y > L/2) for (b) and over 
the lower side (Y < L/2) for (c). 



To have volume-mode bands with finite Chern inte- 
ger, we generally need multiple-band degree of freedom 
within a unit cell. To this end, we considered two peri- 
odic arrays of ferromagnetic particles; decorated square- 
lattice model and honeycomb-lattice model. For the dec- 
orated square-lattice model, we observed that, on in- 
creasing the out-of-plane field, there appears a sequence 
of band touchings between pairs of neighboring volume- 
mode bands. Owing to these band touchings, the Chern 
integers for the volume modes change their signs and, 
concomitantly, the chiral edge mode changes its sense of 
rotation from clockwise to anticlockwise or vice versa. 
For the decorated honeycomb-lattice model, we observed 
a finite band gap between the lowest spin-wave volume- 
mode band and second lowest spin- wave band which 
are connected by a chiral dispersion of an edge mode. 
Though its sense of rotation being unchanged by the 
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FIG. 16: (Color online) Spatial-resolved Fourier power spec- 
tra of magnetization dynamics in the presence of stronger dis- 
sipation. Spatial distribution of the intensity at a; = 29GHz, 
which is obtained by the application of the pulse field at the 
edge, (a) a = 0.001 (b) a = 0.01. 



strength of the field in the honeycomb lattice case, the 
gap and the chiral edge mode persists for a sufficiently 
large field. 

To interpret these results, we next construct tight- 
binding descriptions for the linearized Landau-Lifshitz 
equation, in which atomic orbitals such as 5- wave, p±- 
wave and (ia.2_^2-wave orbitals are introduced within each 
unit cell. Among other, complex- valued characters in the 
p^-wave orbitals break both the time-reversal symmetry 
and mirror symmetries of the models. These symmetry 
breakings lead to a non-zero Chern integer for spin-wave 
volume-mode bands and associated chiral spin- wave edge 
modes. Using this tight-binding model, we argue that the 
level inversions among different atomic orbital levels give 
rise to the so-called inverted spin-wave bands with non- 
zero Chern integers. Our tight-binding analysis for the 
square-lattice model gives quantitative criteria for the 
emergence of finite field ranges within which spin-wave 
volume- mode bands have non-zero Chern integers. 

For the decorated honeycomb lattice model, we em- 
ploy a perturbation analysis, starting from the large field 
limit. The analysis suggests that the effective Hamilto- 
nian in the large field limit always respects time-reversal 
symmetry and the hexagonal symmetry within the order 
of 0(1). Due to the mirror operations in the hexago- 
nal symmetry, the lowest two spin- wave bands form gap- 
less Dirac cone spectra at two inequivalent i^-points. 
Once 0(l/i^)-order corrections are included, however, 
the time-reversal symmetry is lost and hexagonal sym- 
metry reduces to its abelian subgroup having no mirror 
symmetries. As a result, the gapless Dirac-cone spec- 
tra acquire a finite mass of the order of 0{1/H), which 
leads to non-zero Chern integer for the two lowest spin- 
wave bands. This argument explains why the spin- wave 
volume-mode bands with non-zero Chern integers and 
associated chiral spin-wave edge mode persists in a very 
wide range of the field in the decorated honeycomb lattice 



model. 

Since a state-of-the-art sample production does not 
necessarily guarantee perfect periodic structurings, con- 
sidering disorder effects associated with the lattice peri- 
odicity are experimentally relevant, which can be specu- 
lated from well-established knowledges of integer Quan- 
tum Hall physics The effects are two-folded. 
When the strength of the disorders is smaller than a char- 
acteristic frequency scale of the band gap, those volume 
modes near the band gap become localized due to the 
disorders, while chiral edge mode itself is free from these 
weak disorders. As a result, the frequency regime for 
the chiral spin- wave edge mode becomes even wider than 
that in the clean limit. When the strength of the disor- 
ders is increased to be larger than the scale of the band 
gap in the clean limit, however, the 'mobility gap' closes 
and reopens. After the reopening the gap, the topological 
chiral edge mode disappears J^i^iiS^ The proposed chiral 
spin-wave edge mode is also robust against the bound- 
ary shape; the edge modes persist in almost arbitrary 
shapes of the boundary, provided that the edge mode in 
the boundary has no interference with the other mode 
running along the opposite sample boundaryJ^ii^ 

It is also a non-trivial issue whether submicrometer- 
scale ferromagnetic islands behave as a single spin or 
not. In preceding experimental systems mentioned be- 
fore ^i^i non-isotropic shapes of ferromagnetic islands 
give rise to strong magnetic dipolar anisotropics, forcing 
all the spins in each island to point along a same direc- 
tion. In our model calculations, magnetic anisotropics 
within each island are not included from the outset. It is 
interesting to include these magnetic dipolar anisotropics 
into the present Landau-Lifshitz equation phenomenolog- 
ically as the 'single-ion' type magnetic anisotropics. It is 
also equally likely that a ferromagnetic island has a cou- 
ple of low- frequency relevant modes having different spin 
textures within the island. Such modes can be also uti- 
lized as a kind of 'atomic orbital', so that a system with 
only one ferromagnetic island within a unit cell could also 
have a chance to provide a volume-mode bands with fi- 
nite Chern integers. Exploring such systems is, however, 
beyond the scope of the present paper and we leave them 
for future open issues^^ 
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